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Presently, there are several experimental setups dedicated to rare event searches, such as dark 
matter interactions or double beta decay, in the building or commissioning phases. These experi- 
ments often use large mass detectors and have excellent performance in terms of energy resolution, 
low threshold and extremely low backgrounds. In this paper we show that these setups have the 
possibility to exploit coherent scattering on nuclei to detect neutrinos from galactic supernova ex- 
plosions, thus enlarging the number of early detection "observatories" available and helping in the 
collection of valuable data to perform flavour-independent studies of neutrinos' emission spectra. 



I. INTRODUCTION 

Neutrinos from core-collapse supernovae are messen- 
gers of rich information in both particle physics (neu- 
trino properties, osciUations) and astrophysics (super- 
nova mechanism, very dense matter behaviour). They 
also constitute, as of today, the only prompt detectable 
signal of a supernova event, as the technology to detect 
gravitational waves is still under development and no sig- 
nal has been detected yet. 

Charged current (CC) scattering based experiments 
such as Super-Kamiokandc [1], Borexino [2] and LVD [3] 
arc able to detect incoming electron antincutrinos {D^,) in 
the supernova energy range with high efficiency by means 
of the inverse beta decay on free protons (the Chercnkov 
or scintillation light produced by the positron emitted 
during the process is the actually detected signal). Elec- 
tron neutrinos {v^) and (sum of i'^, i/^, Vt, Vr) can 
also be detected by these experiments, but the cross sec- 
tions of the involved processes (CC and neutral current, 
NC, scattering on electrons) are much smaller. Since 
the cross sections for neutrino-electron interactions, es- 
pecially in the MeV energy range, are small, the detec- 
tor mass has to be overwhelmingly large (0(1000 tons)) 
to compensate. Furthermore, the detection capability is 
almost entirely limited to electron neutrinos, while dur- 
ing a supernova explosion neutrinos of all three flavours 
are supposed to be produced. Hence, the inclusive de- 
tection of all neutrino species could provide important 
and oscillation-independent information about their to- 
tal emission flux and spectrum. 

A very promising but not yet exploited mechanism to 
detect Vx is neutrino-nucleus coherent elastic scattering 
on target nuclei. This process is flavour-blind and, for 
small enough momentum transfer, the cross section is 
highly enhanced by the coherent superposition of inter- 
action probabilities for all nucleons within the scattered 
nucleus. Due to the possibility of detecting all neutrino 
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components and the enhancement of involved cross sec- 
tions, the expected number of events from a standard 
supernova turns out to be large enough to make a 1 ton 
scale detector based on coherent scattering as effective as 
a 100 ton light water Chercnkov detector. Moreover, the 
recoil energy of coherently scattered nuclei is correlated 
to the neutrinos' energy in such a way that some infor- 
mation about the neutrinos spectra, the average temper- 
ature for example, can be reconstructed. This is not pos- 
sible, for instance, with inverse beta decay interactions in 
scintillation detectors, where the measured deposited en- 
ergy does not depend on neutrino energy as long as it is 
above threshold. A large mass coherent scattering detec- 
tor can therefore be used, in principle, as a thermometer 
for Vx emitted by collapsing stars. 



Demonstrating the capability of an experiment using 
coherent elastic scattering to detect supernova neutrinos 
increases the number of experiments potentially involved 
in early supernovae detection networks like SNEWS [4]. 



Presently, many experiments for rare events (double 
beta decay, dark matter search) are in building or com- 
missioning phase. These experiments, which are often 
based on cryogenic detectors, have in common good en- 
ergy resolution (hence low threshold capabilities), ex- 
tremely low background and large masses, and they often 
use detectors containing high atomic mass elements (Ge, 
Te, Cd,W). Noble gases (Ar, Xe) and large mass stan- 
dard scintillating detectors (Nal) are interesting as well. 
In Table I some experiments that could potentially use 
this technique are reported. The purpose of this paper is 
to show a systematic study of the potential that different 
materials have as targets for coherent scattering inter- 
actions as a function of the target properties and the 
neutrinos spectra. Materials already used (or planned to 
be used) in large mass rare event detectors are especially 
considered. 
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Experiment 


Detector material Mass [k; 


GERDA (phase II) [5] 


Ge 


37.5"- 


SuperCDMS (phase B)[6] 


Ge 


145"- 


CUORE[7][8] 


TeOa 


74lb 


COBRA[9] 


CdZnTe 


0.42 


CRESST[10] 


CaW04 


10 


XENON100[11] 


Xe 


62'' 


WARP [12] 


Ar 


150 


DAM A/LIBRA [13] 


Nal 


250 



" The experiment feasibility has been demonstrated (project) 
^ The experiment is in commissioning phase 



TABLE I. Experiments that could potentially be able to de- 
tect supernova neutrinos through coherent scattering. 



II. THEORETICAL BACKGROUND 

A. Type II SN 

A core collapse supernova (or type II supernova) is an 
astronomical phenomenon marking the end of a massive 
star's life. Models have shown that for stars vifith masses 
grater than '^9 solar masses the end of the hydrogen 
burning phase is followed by a series of predictable cy- 
cles of contraction, heating and burning of progressively 
heavier elements within the star core (which assumes a 
onion like structure with the heavier element, iron, at the 
centre). The dynamical stability is granted, in each layer, 
by the energy produced in the nuclear fusions. However, 
in the iron core no net energy is produced as no fusion 
can occur; electron degeneracy pressure is the only force 
that prevents the core from collapsing. When the Chan- 
drasckhar limit (1.4 Msun) is exceeded, gravity becomes 
stronger than electron degeneracy and the iron core col- 
lapses. During the collapse, the temperature and density 
increase dramatically and two phenomena occur: 

Phot odis integration 

{A, Z) + J ^ {A' , Z') + {Z - Z')p+ (1) 

{A-A'-Z + Z')n 
Inverse Beta Decay 

+ p — S> n + (2) 

+ n ^ p + De 

The iron core continues to shrink until the density ap- 
proaches the nuclear density, and strong forces and 
neutron degeneracy prevent further collapse. The in- 
falling matter rebounds, creating an outgoing Shockwave 
that dissociates nuclear matter, losing energy and finally 
stalling. The interaction between the Shockwave and the 
core generates extreme temperature/density conditions 
where nucleon bremsstrahlung and pair annihilation take 



place. 

Bremsstraiilung 

{A, Z) + 7 ^ l^e,p,r + i^e,M,r + {A, Z) (3) 

Pair Annihilation 

e+ + -5- iye,^,T + i>e,^i,T (4) 

These last two (Equations 3 and 4) are Zq mediated neu- 
tral current processes. 

Numerical simulations [14] show that the interaction 
of a small fraction (0.1%) of the neutrinos generated in 
this phase with the nuclear matter behind the stalled 
shock should be enough to rise the shock total energy 
to positive values. Unbounded layers are ejected in the 
supernova explosion. 

The processes in Equations 2, 3 and 4 are the mech- 
anisms that generate the neutrino fluxes emitted in the 
supernova explosion. The individual contributions to the 
total flux are '^10-20% from inverse beta decay and ^80- 
90% from pair annihilation and bremsstrahlung [15]. 

A simplified model for the neutrino emission is used in 
literature [16, 17] when the process described in Equa- 
tions 3 and 4 can be considered as the main channels 
of neutrino production. This approximation is especially 
valid in the case of detection through coherent scattering. 
As will be explained in section II B, coherent scattering 
is blind to neutrino flavour; processes like the ones in 
Equation 4 can thus be considered the main source of 
the interacting neutrinos. 

The same simplified model predicts the equipartition 
of the total energy (~ 3 x 10^'^ ergs) among the six neu- 
trino and antincutrino flavours at production. The emis- 
sion spectra will have different shapes due to the dif- 
ferent interaction cross sections, free paths and neutrino 
spheres' radii for the different species, with resulting dif- 
ferent temperatures for i/g, and v^, where are all 
the remaining neutrino and antincutrino (z/^. j/^, i^t) 
families (see Figure 1). Boltzmann spectra with differ- 
ent temperatures are an adequate approximation for our 
purposes and are used in the literature as well [16, 18]. 

The spectra in Figure 1 are calculated for a source lo- 
cated at a distance of 8.5 kpc. This is a common assump- 
tion in the literature [19-21], as it is the distance of the 
centre of our galaxy. Sometimes 10 kpc is used, as it can 
be calculated [22, 23] to be the distance with the highest 
probability of a supernova collapse occurrence. 

B. Coherent scattering on target nuclei 

Coherent nuclear elastic scattering is a neutral current 
weak interaction. From a theoretical point of view, it 
is the same process of neutrino-nucleon neutral current 
scattering. If the momentum of the incoming neutrino 
is small enough, the single nucleon components (protons 
and neutrons) will not be distinguished and the nucleus 
will be scattered as a whole. The scattering amplitudes 
for the different nucleons then coherently sum to give 
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FIG. 1. Boltzmann spectra (at a distance d=8.5kpc from the 
source) for the three neutrino families used in the calculations. 
Green = (3.5 MeV), blue = Ve (5MeV), red = (8MeV) 
and black = total spectrum. 



the total cross section. This turns out to be enhanced 
by a factor of the order of the square of the neutron 
number compared to that of a single nucleon. The co- 
herent behaviour of the interaction will depend on the 
actual momentum transfer between the incoming parti- 
cles (neutrino and nucleus). The higher the momentum 
transfer, the higher the capacity of the neutrino to dis- 
tinguish the single components of the nucleus, and hence 
the smaller the cross section. The nuclear form factor 
(see section IIIB) is the parameter that accounts for this 
dependence of the cross section on the momentum trans- 
fer. 



III. EXPERIMENTAL IMPLICATIONS 

As the aim of this paper is to demonstrate the feasibil- 
ity of neutrino detection through nuclear coherent scat- 
tering, the expected signal for a given detector material, 
mass, neutrino fluxes and spectra have been calculated. 
The number of interacting neutrinos depends on: 

• neutrinos' flux 

• neutrinos' energy spectra 

• target nuclear properties 

• number of target nuclei 

In the following analysis, 1 ton of material is considered 
(TOdctoctor = 1 X 10^ g) and the corresponding number of 
target nuclei is calculated, taking into account stoichio- 
metric ratios of the different atomic species 



(5) 



where a runs over the nuclear species, Aa is the atomic 
mass (the average over the various isotopes is consid- 
ered) , rja is the stoichiometric ratio of the corresponding 
atom and Na is Avogadro's number. 



A. Neutrino properties 

As described in section II A, the largest fraction of the 
energy emitted by a type II supernova (^ 3 x 10^^ ergs) is 
carried by neutrinos. The number of neutrinos depends, 
of course, on the mean energy and the spectral shape. 
In general, the flux of neutrinos of a given flavour i = 
Vp.Up, Vt at a distance d from the source is 



(6) 



where 0^ {E) is the emission spectrum at the source. This 
is assumed to be isotropic. 

Though it is obvious, it is still very important to stress 
the strong dependence of the flux (and consequently 
the expected number of interacting neutrinos) on the 
distance at which the supernova occurs. Calculations 
[22, 23] show that the probability distribution for the 
distance between the Earth and a supernova does not 
exclude smaller distances at all. 

The spectrum (j>i{E) can be considered, as a good ap- 
proximation, to be Boltzmann shaped (Figure 1), 
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where Ti represents the different temperatures {ksT ~ 
3.5,5 and 8 MeV) and Ni is the total number of radiated 
neutrinos. Assuming the equipartition of energy, it is 
N^^ = 3.0 X lO^*^, TVp^ = 2.1 X 10^6 ^nd N^^ = 5.2 x 10^^. 

Importantly, possible effects of neutrino oscillations 
have been neglected in this paper. However, this should 
not introduce any significant error in the signal estima- 
tion as the coherent scattering is a flavour-blind pro- 
cess which could, in principle, lead to an inclusive and 
oscillation-independent detection of all (non-sterile) neu- 
trinos. The spectra of the different neutrino families can 
then be summed in a total spectrum. 



(8) 



2 — l^e I'^e -i^x 



Its shape is not analytically defined but strongly depends 
on the temperatures of the single families. The depen- 
dence of the final result (i.e. the amplitude of the de- 
tectable signal) on the spectral parameters will be ana- 
lyzed in section IV. 



B. Target properties 

The cross section for the interaction between incoming 
neutrinos with energy E and target nuclei via coherent 
elastic scattering is [19, 24] 
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where Gp is the Fermi constant and 6 is the angle be- 
tween the original and the scattering directions. is 
the weak charge of the nucleus. This last factor is the 
one that accounts for the enhancement of the cross sec- 
tion due to the coherent superposition of single-nucleon 
cross sections 
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Q^,=N- {I - Asin^ ew)Z 



(10) 



where N and Z arc respectively the number of neutrons 
and protons within the nucleus and sin^ Qw ~ 0.231, 
which means that almost only neutrons contribute to the 
weak charge. The last term in Equation 9 is F{Q^), and 
it is the elastic form factor at momentum transfer Q 



Q2 = 2£;2(1 -cos6i) 



(11) 



and represents the distribution of the weak charge within 
the nucleus. The proton density distribution is often well 
constrained by measured charge densities, and models 
exist to calculate the overall form factor. Probably the 
most complete treatment is in [25-27] and the form factor 
used is 



(QRo) 



X exp 



QRo 



QRq 



(12) 



where Rq is the nuclear radius defined as 



Rl = R^ - hs^ 
R = (1.2 X Ai/3)fm 

and s the nuclear skin thickness (0.5 fm). Figure 2 shows 
the form factor for some of the nuclei considered in this 
paper. The heavier the nucleus, the more important 
the correction introduced: the incoming neutrino will no 
longer coherently sec the nucleus as a whole at a smaller 
momentum transfer if the nucleus is larger, while for 
smaller nuclei the coherent behaviour lasts until larger 
values of the momentum transferred in the interaction 
(corresponding to a larger recoil energy of the nucleus) . 



C. Events number calculation 

Inserting Equations 10 and 12 into Equation 9, the dif- 
ferential cross section for the cited nuclei can be com- 
puted for an interacting neutrino of a generic energy 
(as an example. Figure 3 represents cross sections for a 
50MeV neutrino). 

To obtain the number of interactions within a certain 
mass of material, the cross section times the neutrino flux 
has to be integrated. As the cross section depends on the 
scattering angle, and the scattering angle determines the 
recoil energy of the scattered nucleus {2MT = ~ 
2E^{1 — cos 9), M is the mass of the recoiling nucleus) 
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FIG. 2. Form factors for different nuclei. 
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FIG. 3. The polar plot of the coherent scattering differential 
cross section for different nuclei for a 50 MeV neutrino shows 
the angular dependence. The colour coding is the same of 
Figure 2. 



through the momentum transfer Q, the events' yield can 
be obtained through a numerical integration as a function 
of the kinetic energy of the recoiling nucleus T. The 
analytical form of the integral is 



Y{T) 



diV.„ 



dT 

= E 

a— nuclei 



dQdEx 



X 



-^(Q\F{Q'),Qi^,A)^UE)x 
n2 



dil 
xs(t 



2Ma 



(13) 



where is the total number of target nuclei 

(Equation 5) and the sum runs over the different nuclear 
species in the detector material. The three different neu- 



5 



trino species are already summed using (f)^^^ . In the case 
the detector contains a single nuclear type (as is the case 
in Ge or noble gas detectors), the first sum is redun- 
dant, while for compound materials (Te02, Nal) the to- 
tal number of target nuclei for each nucleus depends on 
the stoichiometric ratio. 

The integration has been performed numerically for 
each neutrino type independently and for the differ- 
ent nuclei. In Figure 4 the yield for 1 ton of BGO 
(Bi4Ge3 0i2) scintillator is presented as an example. 




Recoil Energy [keVJ 



FIG. 4. Events yield for 1-ton of BGO scintillator experiment. 



The different behavior of light and heavy nuclei is evi- 
dent: the bismuth cross section is enhanced by the large 
number of neutrons leading to a high number of inter- 
actions, but the energy of the recoiling nucleus is always 
small due to the large mass. Hence, the spectrum of 
nuclear recoils is steeply exponential shaped. Moreover, 
high energy recoils are further reduced by the form factor 
that damps the cross section at high momentum transfer. 
On the contrary, the number of oxygen recoiling nuclei 
is much smaller but almost constant as a function of the 
recoiling energy. 

The obvious consequence of the recoil spectrum shape 
is the importance of the energy threshold and perfor- 
mance of the detector close to it, both in terms of ef- 
ficiency and background. This is the reason why rare 
events experiments (that usually have low threshold ca- 
pabilities and very low and well known background), in 
particular dark matter and double beta decay bolometric 
experiments, have the potential for exploiting coherent 
scattering as a flavour-blind supernova neutrino detec- 
tion mechanism. 

In Figure 5, the response function (differential events 
yield. Equation 13, times the recoil energy) for 1 ton of 
material is depicted, while in Figure 6 the total number of 
events above threshold as a function of energy threshold 
is calculated for the same compounds. 
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FIG. 5. Events yield for 1-ton of Nal, Te02, Xe, Ge, CaW04 
and BGO. 
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FIG. 6. Number of events above threshold as a function of 
the energy threshold. 



IV. UNCERTAINTIES PROPAGATION 

The main uncertainties in the results reported in 
section III are due to the propagation of uncertainties in 
the astrophysical models of supernova explosions, hence 
in the emission spectra of neutrinos. If the equipartition 
of energy among the species is a widely accepted state- 
ment deriving from the universality of the interactions 
involved in Equation 3 and Equation 4, parameters such 
as the total energy of the explosion and the average en- 
ergy, or the temperature, of the neutrinos spectra are 
much more uncertain. 

A variation in the total energy carried by the neutrinos 
(assuming a simple luminosity scaling and no change of 
spectral shape) has the trivial effect of changing linearly 
the number of neutrinos and hence the number of inter- 
actions in a given detector. The effect of a variation of 
the temperature of the spectra described in Equation 7 
must be discussed in more detail. A lower temperature 
corresponds to a red-shifted emission spectrum. Since 
the total energy content is the same, but the average en- 
ergy is smaller, the number of emitted neutrinos is larger 
but a larger fraction of them will not be able to produce 
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recoiling nuclei above threshold energy. Hence, the signal 
in the detector drops significantly. 



30- 
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FIG. 7. Number of events in 1 ton of BGO scintillator with 
3keV threshold for different values of neutrinos spectra tem- 
perature parameter. 



An increase of the temperature, on the contrary, leads 
to a higher maximum transferred energy between the in- 
coming neutrinos and the recoiling nuclei. However, since 
the total number of neutrinos is smaller the signal in- 
crease is dumped and tends to saturate. 



V. CONCLUSIONS 

Reported calculations show that a 1 ton detector with 
adequate energy threshold (as low as a few keV) has the 
potential to detect supernova neutrinos through coher- 
ent scattering on nuclei. The sensitivity of such an ex- 
periment depends on another very important parameter, 
the background. The number of signals coming from in- 
teresting events has to be compared to the number of 
events coming from the background in order to deter- 
mine the sensitivity. Rare events experiments usually 
have very low backgrounds and many of the detectors 
used in this field are able to distinguish heavy particle 
(nuclear) recoils from electron recoils through simultane- 
ous measurement of different energy channels, often be- 
coming zero-background experiments. As recently sug- 
gested in [28], these experiments' capabilities could not 
be limited to supernovae neutrino detection, but extend 
to low energy neutrino measurements in a wider frame- 
work including short baseline oscillation studies with very 
intense sources. 



The resulting propagation of the uncertainties on the 
temperature to the signal in the detector is reported in 
Figure 7. A very large variation of the Vx temperature 
is considered while the ratio between the temperatures 
of the three different species of neutrinos has been kept 
constant. The result is the number of events in 1 ton of 
BGO scintillator with an energy threshold of 3 keV. The 
shadowed region represents a reasonable uncertainty of 
30% on the temperature parameter. 
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